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Abstract 



It is shown that that violation of causality in two-dimensional light- 
front field theories quantized in a finite "volume" L with periodic or 
antiperiodic boundary conditions is marginal and vanishes smoothly 
in the continuum limit. For this purpose, we derive an exact integral 
representation for the complete infinite series expansion of the two- 
point functions of a free massive scalar and fermi field for an arbitrary 
finite value of L and show that in the L — > oo limit we retrieve the 
correct continuum results. 
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1 Introduction 



Light-front field theory has some important advantages ^[^HH] due to it sim- 
plified vacuum structure that make it a promising theoretical framework for 
elementary particle physics. On the other hand, a systematic understanding 
of non-perturbative structure of light-front (LF) field theory has so far not 
been achieved. A particularly advantageous formulation appears to be LF 
quantization in a finite-volume, known as discretized light-cone quantiza- 
tion (DLCQ) |U[3]. It incorporates in an efficient way boundary conditions 
which are required for mathematical consistency even in the continuum for- 
mulation The point is that the numerous LF constraints which reduce 
the number of independent field degrees of freedom can be (at least in prin- 
ciple) uniquely inverted only if the corresponding Green's functions satisfy 
(anti)periodic boundary conditions. For periodic boundary conditions, one 
may then study physical implications of zero-mode operators (carrying the 
LF momentum p + = 0) in a finite volume which serves as infrared regular- 
ization. Thus, it is desirable to formulate light-front quantization in a finite 
space volume with a discrete infinity of modes as a quantum field theory in 
its own right. This implies that one should verify all usual well-established 
properties (such as causality, Poincare invariance, singularity structure, etc.) 
in this framework to check its overall consistency. 

It is far from clear that these properties will hold true. Indeed, already 
for the simple question of microcausality of a massive scalar field in two 
dimensions it has been concluded that causality is violated by the infrared 
finite- volume regularization [Jj. More precisely, it has been argued that 
periodic boundary conditions are incompatible with causal behavior of the 
light-front quantum theory. The method used to demonstrate this was basi- 
cally a numerical study of the corresponding series truncated at some value of 
discretized LF momentum p + for which the results stabilized. This method 
gave a very satisfactory picture of the causality in a space-like box, namely 
vanishing (up to negligible numerical effects) of the Pauli-Jordan (PJ) func- 
tion (which is twice the imaginary part of the full Wightman function) for 
space-like separations and a usual oscillatory behaviour in the time-like re- 
gion. This picture however failed completely for a LF system restricted to 
— L < x~ < L. Not only did the numerical results for the PJ function fail 
to vanish for x 2 < 0, but it was even found not to converge to the correct 
continuum expression. Two obvious explanations are: 1. the discretized 
light-front theory has some fundamental difficulty or, 2. the results of the 
numerical computations are misleading or at least do not reveal the full na- 
ture of the problem. To clarify this issue it would be preferable to find a 
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method for analytical evaluation of the infinite series expansion of the PJ 
function, corresponding to the integral representation of the PJ function in 
the continuum formulation. 

A calculation of this kind has been sketched in Ref.[S] indicating that 
the PJ function computed in a finite volume converges to the correct con- 
tinuum form in the large L limit. The results of Ref. 7 , however, put this 
picture into doubt (see also Ref. 8^ for a careful numerical analysis). We 
find it very important to clarify the situation by an independent analytical 
calculation and this is the main purpose of the present paper. More gen- 
erally, we wish to show that there is a natural mathematical formalism for 
evaluation of infinite series corresponding to various correlation functions 
of the discretized LF theory. This formalism is very well adapted to the 
form of LF kinematics and dispersion relation and it uses some properties of 
polylogarithm functions ^U]. As a result, an integral representation can be 
given for the infinite series expansion, in particular for that corresponding 
to the free-field Wightman functions. This representation explicitly depends 
on the box length L. We use analytical methods to study the large L de- 
pendence of the integral representation of the PJ function. The result of 
our analysis is unambiguous: there is no physically relevant violation of 
causality. We recover the well known continuum result plus finite-volume 
corrections which are suppressed by a power of L and thus vanish in the 
large- L limit. Moreover, even for relatively small value of L one obtains a 
very plausible picture: an oscillatory structure in the time-like region and 
near-vanishing values of the PJ function in the space-like region. This will 
be demonstrated numerically in Sec. 4 below. 

2 Free-field correlation functions 

We first briefly describe the derivation of the two-point Wightman function 
for a massive LF scalar field in two dimensions. In the continuum formula- 
tion, the mode expansion of the scalar field is 1 

oo 

<f>{x) = J ^L-[ a ( q + }X +) e -h + *- +o t (g + ,x + )e3« +x ~]. (1) 



1 We use the convention x — (x + ,x~),x 2 — x + x~,d^ = 2d± = d/dx ± 
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In the free case, the time evolution is simple. From the Klein-Gordon equa- 
tion we have a(q + ,x + ) = a(q + , 0) exp(— ip~x + /2), p~ = fi 2 /p + , and hence 



4>(x) 



dq+ 



[a(q + )e 



(2) 



where a((7 + ) = a(q + ,0) and small imaginary parts for x are understood 
(see below). The corresponding two-point function 



D{x) = <O|0(O)^(x)|O) 



dp + 
4irp^ 



(3) 



(4) 



can easily be obtained using the commutation relation 
[a(p + ),a\p' + ] = 4irp + 5(p + -p' + ). 
With the help of known integral formulae 3.471 from one further finds 

D(x) = -^|e(x+) + e(x-)|(^ (/wv / ^) -i e(x + )J (/iV / ^2)) 

(5) 



+ "He 2° 

47T 



(x-)\K (fiV^), 



where e(x ± ) is the sign function. To guarantee the convergence in equations 
01 and El below, it is assumed that the quantities x^ include an infinitesi- 
mal positive imaginary part. Jo(z), Kq(z), and Nq(z) are standard Bessel, 
modified Bessel, and Neumann functions. 

The corresponding expression for the same system quantized in a finite 
volume with field obeying periodic boundary condition reads 



D(x) 

where the field expansion 



1 ^ — ^ 1 2^ >n ' C 2 — T" 3 '"'" -t- 

4vr H n 

n=l 



2vr 

T 



-n. 



(6) 



4>(x) 



l ^ l 



E 



n=l yp n 



with 



2 2 

_i r ,+ r -_ L^ T + i n + r --i_i ii^ T + 

a n e p ™ + w n e p ™ 



(7) 



(8) 
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has been used. The mode with n = is excluded in the above series since 
the free Klein-Gordon equation (4<9 + <9_ + /i 2 )0(x) = is not compatible 
with a non- vanishing Fourier mode carrying p + = if ji ^ 0. 

A similar treatment can be given for free massive LF fermions. In the 
representation where the 75-matrix is diagonal, the LF Dirac equation for 
the spinor field ip T = (tpi, ip 2 ) of mass m separates into two equations for 
the components ipi,ip 2 : 

2id + ip 2 (x) = mipi(x), 2id-ipi(x) = mip 2 (x), 75 = diag(-l, 1) (9) 

The independent component V'2 can be expanded at x + = into the Fourier 
integral with the operator coefficients b(p + ), b^(p + ), d(p + ), cft(p + ) and the 
above dynamical equation determines its LF time evolution as 

00 , + . 2 .2 

Mx)= I P ^ [b(p + )e~^ P+x ~~^ x+ +d\p + )e^ P+x ~ + ^ x+ }. (10) 
J A7r\/p+ 


The solution of the constraint equation is 

°° A + 2 2 

Mx) = m [ — ^^—[b{p + )e~ ip+x ~~^ x+ -d\p+)e hp+x ~ H ^ x+ }. 

J 47TVP+P+ 

v 

(11) 

The quantization rule 

{V> 2 (0, x~),4(0, y-)} = l -5(x- - y-) (12) 

is equivalent to the following anticommutation relations for the Fock oper- 
ators: 

{b(p + ),tf(p' + )} = {d(p + ),d\p' + )} = ^p + 5{p + -p' + ). (13) 
It is very simple to calculate the two-point Wightman functions S a p{x): 

S 22 ( x ) = (0|V 2 (0)4(x)|0) = J ^-e ip+x ~ + ^ x+ 



Su(x) = (0|Vi(0)^(x)|0) = m 2 f JV e ^ + *- +1 2^ + 

J 47TO+ 


°° A + ■ 2 

S 12 (x) = (0|Vi(0)4(x)|0) =mj ^e^ p+x ~ + ^ x+ . (14) 
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Using again the integral formulae we obtain the explicit expressions 



S22(x) 



Su(x) 



Su(x) 



|e(x + ) + e(x~ 
(e(x + ) — e(x 
|e(x + ) + e(x~ 



,m x 




8 V x 



m x tz i I 2^ 

— \\ — Ai mv — x z ) 

4tt V x 



Ji(mV x 2 ) + i e(x )Ni(mVx 2 ) 



,m x 




J\(mv x 2 ) + i e(x + )Ni(mv x 2 ) 



+ {e{x + ) — e(x )j -j—v —Ki(m\/—x 2 ) 



e(x + ) + e(x~)|— No(mV x 2 ) — i e(x + )Jo(mV x 2 ) 



+ |e(x + ) — e(x )\^Ko(m\/—x 2 ) 
The corresponding expansions in the discretized case are 



(15) 



tpi(x) = - 



n=l 
m \ - 1 



V — 

'2L Pn 



b„e 



.iv+ x - -±™ 

2 Pn 2 7T 



9 I J n o + 



+ die 2 



die 2 



where the fermion Fock operators satisfy 



(16) 



(17) 



To keep our discussion close to the scalar field case, we will work with 
periodic fermi field ip(—L) = ip(L) (it is straightforward however to perform 
the whole analysis for antiperiodic fermi field). The field equations © again 
require zero modes of both components to vanish. 

The above field expansions yield for the correlation functions 



522(2;) 
Sii(x) 

£12(20 



2L ^ 

n=l 

m 2 L ^2, 1 y+x-+i?4 



E 



8vr 2 — n* 



n=l 

00 



1 hvlx 



4tt , n 

n=l 



— e 



(18) 
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Do the discrete representations (jHJ) and (|18|) lead to the continuum results 
© and ()15|) for L —* oo? More specifically, does 2 Iml) reproduce the 
continuum Pauli-Jordan function in this limit? To answer this question, 
we begin by replacing the infinite series © and ()18|) by integrals using an 
integral representation of the polylogarithm functions. 

3 Integral representation of the discrete two-point 
function and the causality problem 

We start with the discretized two-point Wightman function It is ex- 
pressed as an infinite series. A very useful alternate representation of this 
function can be given in the form of an integral based on polylogarithms. 
Consider the function of two independent complex variables defined by 

oo n 

Fi(z,q) = e«/». ( 19 ) 

For any finite q it can be shown that the power series converges only within 
the unit circle \z\ < 1. Expanding e q l n in powers of its argument we obtain 

*M = E { f- E = E fr (20) 

A;=0 n=l fe=0 

Here we ticivG used, th.6 definition jlQj of the polylogarithm function Li m , 

oo n 

= E ^r- ( 21 ) 

n=l 

Note that this series representation ()21j) of Li m converges only if |z| < 1. 
Its analytic continuation to the rest of the complex z plane is provided by 
the integral representation 

oo 

Lir n (z) = J ^jdu^- l ,(rn>l), (22) 

which shows that Li m (z) is actually analytic throughout the z-plane except 
for a cut on the positive real z axis, linking +1 to +oo. Substituting this 
formula into (|2Uj) and interchanging integration on u with the summation 
on k, we arrive at the integral representation of the series (|19|) . 

roc i 

F 1 (z,q)= du M 2y/qu). (23) 

Jo z L e u — 1 
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To obtain this result we have used the identity 



[v) k 



k=0 

Comparing © and (fTU|) we note that 



E^ = 'o(2v^). (24) 



D(x) = ^-F 1 (e i ^,iQ), (25) 

4-7T 



where 



£ = fi 2 x + x-/4, Q=^± x + . (26) 

47T 

This is the first step of our analysis whereby the infinite series © has been 
rewritten in integral form via l)23j) and 1)25 Jl . 

There are four distinct cases to consider each associated with a quadrant 
of the x + ,x~ plane. Consider first the case where both x^ are positive 
(Q > 0,£ > 0). We may rewrite Eqs.^ and (J2SJ) as 

1 POO 1 

t){x + > 0,x~ > 0) = — - / du 7 r I (2e in / A ^u~). (27) 

4irQ Jo exp(Mu-i£))-l 



The continuum limit is then obtained by considering the limiting value of 
the RHS of Eq.(|27|) for Q — > +oo. We now show that the result is the 
reduced version of Eq.© appropriate to the regime x^ > 0. Before pro- 
ceeding we note that I (2e i7r / 4 ^/u) is an analytic function of u throughout 
the entire finite complex plane since its Taylor series expansion in powers of 
the argument has an infinite radius of convergence. The remaining factor 
in the integrand of 1)27)1 is analytic throughout the u plane with the excep- 
tion of simple poles at the discrete set of points u n = i(2nirQ + £) on the 
imaginary axis, where n is any integer. In view of these analytic properties 
of the integrand we can alter the integration contour without changing the 
value of the integral in (|27|) as long as we avoid crossing through any of the 
singular points u n and maintain the given endpoints. The simplest choice 
of a preferable contour is presented in Fig.l. It consists of the straight- 
line segments (Ci) on the imaginary axis, the semi-circle (C2) centered on 
the pole uq, the line (C3) parallel to the positive real axis, and finally the 
straight-line segment (C4) parallel to the imaginary axis. Note that on C4 
we have Re u = R and we will require that R/Q > aQ where a > 2. 

On the contour C\ we have u = iv, where v is positive real. One imme- 
diately finds that this contribution to D(x) may be written as 

J dvJ (2^)cot((v-0/(2Q))-^Q J dvJ (2y/v), (28) 




S 



where P denotes principal value. Each of the two integrals is real and finite. 
The second term can be evaluated in closed form and the result is 



' v;Tf ^" e -J 1 (2 v ^QTi), (29) 



8vr Q 



and in particular it vanishes in the large Q limit. 

In order to obtain the L — > oo limiting behavior of the first term of (|28[). 
one can replace the cotangent function by the inverse of its argument, so 
that after the change of variable v = w 2 /4 one obtains 



oo _ oo 



i r = j_ r ^ = i^^^. (30) 

4vr ' v-£ 2-k J w 2 -4£ 4 Ul VM y ' 

o o 

On the infinitesimal semicircle C2 defined by u — i£ = ee l9 , —tt/2 < 9 < vr/2, 
with e — > + , we may use the approximation 

7 " \ « (31) 



Replacing further the function Jo by its value at uo and using the relation 
[TT] Jo(2«\/?) = ^o(2\/D) we find that this contribution to D is equal to 
j^Jo(2^). Note that this quantity is independent of L and thus survives 
in the large-Q limit. 

On the horizontal semi-infinite line C3 we may write u = i£ + Q{m + v), 
where v is real, positive. The contribution to D(x) is given by 

aQ 

~h\ dv ~r^[M2VZ + KQ-iQv). (32) 


The asymptotic analysis of Ea. (|32|) in the limit Q — > +00 is very lengthy 
and involves substitution of an integral representation of Jo combined with 
the method of steepest descent. The details of this calculation will be given 
elsewhere [12] • The final result is that the leading behavior of the expression 
1|32J) is given by a term proportional to exp(iQ). This vanishes in the limit 
L — > 00 as long as for any finite L the quantity x + includes a small positive 
imaginary part such that L x Im(x + ) — * 00. This requirement is satisfied 
for example by the choice Im(x + ) = 0(L~ 1 / 2 ). 

Finally, points on the line segment C4 are described by u = R+iv, where 
v is real, positive. The contribution to D(x) from C4 is then found to be 
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dominated by exp — (a — V2a)Q and thus vanishes in the large-Q limit 
since we choose a > 2. 

Combining these results we have 

lim D(x + > Q,x~ > 0) = -\(n (hVo?) - iJoifi^/x 1 )), (33) 

in agreement with Eq.Q. 

It is easy to see from the formula (|23|) that results for finite L for the case 
x + < 0, x~ < (£ > 0, Q < 0) can be obtained from those for x + > 0, x~ > 
by complex conjugation: 

^(e-WI.-ilQl) = [^(e^iQ)]*. (34) 
Likewise, we have for x + > 0, x~ < (Q > 0, £ < 0) 

^) = ^i(e~^/ Q ,zQ), (35) 

and this can be used to generate the results also for the regime x + < 0, x~ > 
(£ < 0, Q < 0) according to 

£(z) = ^i(e^l/l«l,-t|Q|) = [^(e-W^Q)]*. (36) 

The evaluation of F\ in Eq. (|35|) for large L proceeds as above using the same 
multi-component contour except that the semicircle C2 is not applicable 
since the pole uq is situated at The final result is 



lim £(*+ > 0, x- < 0) = J- / = -2-^(2^), (37) 



in agreement with the continuum formula In particular, this result 

means that since the imaginary part of the function D{x) vanishes for space- 
like x 2 , the causality is restored in the infinite- volume limit. For large L the 
leading L-dependent terms are of the order L _1//4 . 

Finally we briefly discuss the correlation functions S22 and Su of l|18jl . 
The first of these is a special case of the function 

OO 

F (z,q) = £ = j^+ q jdu hS0L t (38) 

^ 1 — z J z l e n — 1 x /qu 

n=i n v * 
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while Sn is a special case of 



00 n f 

F 2 {z, q ) = Y.- 2 e q/n = du x U u 

^ n z J z 1 e u — 1 

n=l 



u Ii(2y/qu) 



qu 



(39) 



In the regime x^ > we find that 



lim S 22 (x) = Ae^pfdw^^-^U^) 
Q^+oo J w z — 4t V6 



= —^\—(ji(m\G?)+iNi(m\G?)). 

4 y x \ / 

Likewise, in the same regime, we have 



(40) 



lim Sn (x) = \ / — - ( Ji (mVs^) + iJVi (mV? 

q^ +00 v/ m V x + V 1V ; 1V 



(41) 



These results are in agreement with the continuum formulas 1)15(1. 



4 Numerical results 

In principle, one may try to evaluate the integral ()23|) representing the Pauli- 
Jordan function in finite volume numerically for increasing values of L to 
examine the rate of convergence towards the continuum result. In practice, 
this is rather difficult to achieve since the integrand of the representation 
Q23JI oscillates rapidly due to the presence of the Bessel function Iq . Already 
for relatively small values of Q the amplitudes of these oscillations are so 
large and the spacings of successive zeros are so small that it is very difficult 
to reliably evaluate the integral by standard numerical routines. Since it is 
not our goal to perform extensive numerical analyses in the present work, we 
have computed the integral for a few relatively small values of Q{= 4, 15, 18) 
using an integration method based on Chebyshev polynomials as well as by 
a Clenshaw-Curtis method. For definiteness we set /j?x + = 1 so that the 
corresponding box lengths given by L = 4ttQ are approximately L = 50, 188 
and 226. The results are displayed in Fig. 2. An essential difference in the 
behavior of the Pauli-Jordan commutator function between the space-like 
region (negative values of x~/L) and time-like region (positive x~/L) is 
obvious already for the smallest value Q = 4. It is also evident that for 
larger Q the oscillatory behavior of the continuum curve in the time-like 
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region is resolved with increased accuracy. This is particularly true in the 
interval < x~ / L < 0.5 but the number and position of oscillations is 
semiquantitatively reproduced also for x^/L > 0.5. Although the Pauli- 
Jordan function for finite volumes is not zero in the space-like region, it 
is reasonably close to it. We recall that for our choice of Q values we are 
still very far from the infinite-volume limit so the obtained behavior of the 
Pauli-Jordan function is very plausible and consistent with our analytical 
results. 

5 Discussion 

The Discretized Light Cone Quantization method has led to many inter- 
esting results over the last fifteen years. It is however essential that the 
method satisfy all requirements and principles of a consistent relativistic 
theory. It is rather clear that restricting a system to a finite spatial volume 
and imposing (anti)periodic boundary conditions on quantum fields will gen- 
erally lead to a theory which does not have exactly the same properties as 
the continuum (infinite-volume) theory. Some finite-volume effects may be 
present |13j . However, the crucial point is that the continuum limit of a 
finite-volume theory should recover all necessary properties including e.g. 
Poincare symmetry. The principle of microcausality (or locality) is one of 
the fundamental properties of relativistic dynamics and the DLCQ method 
would face a serious difficulty if it would be in conflict with this principle. 
We have demonstrated in the present work analytically as well as numeri- 
cally that this is not the case. With the infinite number of field modes the 
violation of microcausality in a LF finite volume with periodic scalar field is 
only a marginal effect and continuum results including the causal behavior 
are restored in the L — > oo limit. In practice, the DLCQ calculations of 
mass spectra and wavefunctions are always performed for finite L and with 
a finite number of Fourier modes. At this step, the causality may seem to 
be violated [Jj (see also [H] for a treatment that averages over some range 
of L values and restores the causal behavior in a finite volume). However, 
as physical quantities calculated with the DLCQ method have to be extrap- 
olated to the continuum limit, there is no inconsistency, since, as we have 
shown, the causality is restored there. 

Our conclusion has been obtained by means of a well defined mathe- 
matical treatment leading to an integral representation of the infinite series 
representing the two-point functions for a finite volume. We have extracted 
the -L-independent part of the integral as well as the leading large- L cor- 
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rections. Our results for the PJ function 2ImD(x) are consistent with the 
results of Ref.fSJ. The method used enabled us however to calculate the 
complete Wightman functions of two-dimensional free massive bosons and 
fermions quantized in a finite volume. A detailed discussion of our mathe- 
matical treatment will be published separately |12| . 

Ames Laboratory is operated for the United States Department of En- 
ergy by Iowa State University under Contract No. W-7405-Eng-82. One 
of the authors (L.M.) was partially supported by the Slovak Grant Agency, 
VEGA Grant No. 2/3106/2003. 

6 Acknowledgements 

The authors would like to thank J. P. Vary for support and Peter Markos 
for assistence with the numerical calculations. 

References 

[1] P.A.M. Dirac, Rev. Mod. Phys. 21 (1949) 392. 

[2] H. Leutwyler, J. R. Klauder and L. Streit, Nuovo Cim. 66A (1970) 536. 

[3] F. Rohrlich, Acta Phys. Austriaca 32 (1970) 87. 

[4] T. Maskawa and K. Yamawaki, Prog. Theor. Phys. 56 (1976) 270. 

[5] H.C. Pauli and S.J. Brodsky, Phys. Rev. D 32 (1985) 1993. 

[6] P. Steinhardt, Ann. Phys. 128 (1980) 425. 

[7] Th. Heinzl, H. Kroger and N. Scheu, hep-th/9908173 

[8] D. Chakrabarti, A. Mukherjee, R. Kundu and A. Harindranath, Phys. 
Lett. B 480 (2000) 409. 

[9] S. Salmons, P. Grange and E. Werner, Phys. Rev. D 60 (1999) 067701. 

[10] L. Lewin, Polylogarithms and Associated Functions, North Holland, 
New York, 1981. 

[11] I. S. Gradhsteyn and I. M. Ryzhik, Tables of Integrals, Series and Prod- 
ucts, Alan Jeffrey, Editor, Fifth Ed., Academic Press, New York, 1994. 

[12] L. Martinovic and M. Luban, manuscript in preparation. 



13 



[13] R. Haag, Infinite systems, in: Brandeis Lectures (1971). 



14 



i(7cQ + Q 



Im u 




Re u 



Figure 1: Integration contour chosen for our evaluation of the integral |2 
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Pauli-Jordan function, Q=4 
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Figure 2: The commutator function ImD(x + ,x~) evaluated numerically is 
plotted for three values of Q as a function of s = x~ /L and compared with 
the continuum function 1/4 Jo (s/Ls) (solid line) in the box of the length L = 
AttQ. For simplicity, we chose fi 2 x + = 1. The time- like region corresponds 
to positive values of s, the space-like region to negative ones. 
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